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We present an algorithm for calculation of the Gaussian classical capacity of a quantum bosonic 
memory channel with additive Gaussian noise. The algorithm, restricted to Gaussian input states, 
is applicable to all channels with noise correlations obeying certain conditions and works in the full 
input energy domain, beyond previous treatments of this problem. As an illustration, we study 
the optimal input states and capacity of a quantum memory channel with Gauss-Markov noise [J. 
Schafer, Phys. Rev. A 80, 062313 (2009)]. We evaluate the enhancement of the transmission 
rate when using these optimal entangled input states by comparison with a product coherent-state 
encoding and find out that such a simple coherent-state encoding achieves not less than 90% of the 
capacity. 



I. INTRODUCTION 

A central problem of information theory is to derive 
the capacity of communication channels, which is the 
maximal information transmission rate for a given avail- 
able energy. For quantum channels, which may trans- 
mit quantum or classical information, one can define the 
quantum or classical capacity. The present paper is fo- 
cused on the latter. Just as in classical information the- 
ory, the question of whether or not the classical capacity 
is additive is of central importance. If the capacity is su- 
peradditive, then the asymptotic transmission rate can 
be higher than the maximal rate achievable for a single 
use of the channel. Although some quantum channels 
were proven to be additive [TH3], a counterexample was 
exhibited, proving that this does not hold in general 4J. 
The situation becomes rather different for channels with 
memory, as then the subsequent uses of the channel get 
linked. In particular, if the memory is modeled by noise 
correlations between the uses of the channel, it has been 
shown that the transmission rate may be enhanced by 
using entangled input states [51I19j. 

A particular interest has been devoted to Gaussian 
bosonic channels as they model common physical links 
such as the optical transmission via free space or opti- 
cal fibers. An overview of these bosonic channels can be 
found in |20| . Among them the most studied are the addi- 
tive noise and lossy Gaussian channels (see e.g. [3ij211 22' 
and Refs. therein). For these channels, it was also shown 
that correlated noise may lead to superadditivity as the 
transmission rate may be enhanced by input states with 
some degree of entanglement [T2HI9] . We remark that 
all these studies were restricted to the set of Gaussian 
input states, so that the derived quantity can be viewed 
as a "Gaussian capacity". However, if the conjecture that 
Gaussian states minimize the output entropy of Gaussian 
channels could be proven, then this quantity would turn 
into the actual capacity. For a particular case, i.e., the 
memory less lossy Gaussian channel with vacuum noise, 
this conjecture was proven to hold [3]. In the present pa- 
per, we also restrict ourselves to the set of Gaussian input 



states, that is, we investigate the Gaussian capacity. 

In a recent work, we have evaluated the Gaussian 
capacity of a Gaussian quantum channel with additive 
Markov correlated noise, restricted to a certain input en- 
ergy domain [IB] . Correlated Gaussian noise appears, for 
instance, in the models of downlink communications be- 
tween satellites and terrestrial stations |23j and a simple 
description of such correlations can be provided by an 
underlying Markov process. We found a quantum water- 
filling solution to the Gaussian capacity, similar to the 
classical water-filling that appears when considering par- 
allel classical Gaussian channels [24]. This similarity is 
very surprising, as we must take into account that a part 
of the input energy, in addition to being spent on classi- 
cal modulation, is spent on the preparation of quantum 
information carriers (e.g., squeezed states), something 
which has no classical analog. The notion of quantum 
water-filling appeared for the first time in the discussion 
of the capacity of a memoryless phase-dependent Gaus- 
sian channel |21| . although there the quantum informa- 
tion carriers were considered part of the channel and only 
the energy cost of classical modulation was considered, 
thus making this solution a straightforward analog to the 
classical one. 

In this paper, we present a solution to the Gaussian ca- 
pacity of the Gaussian additive noise channel in the full 
input energy domain, where noise correlations are given 
by stationary (shift-invariant) Gauss processes. We show 
that the method is applicable even to a larger class of 
noises. We present an algorithm for a numerical solution 
of the arising optimization problem based on the method 
of Lagrange multipliers. Although the algorithm was de- 
rived independently, the validity of this method is based 
on arguments which are essentially equivalent to those 
presented in [T7] for a lossy Gaussian channel where the 
quantum water-filling solution was also obtained. Us- 
ing this method, we analyze the Gaussian capacity and 
the associated optimal input states and encoding as a 
function of the noise parameters for the special case of 
Gauss-Markov noise, including the limiting case of maxi- 
mal noise correlations. In addition, we evaluate the gain 
from using the optimal input states, which are entan- 
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gled and therefore may be complicated to produce, with 
respect to easily generated coherent product states. 

The paper is organized as follows. First, we introduce 
the notion of classical capacity and Gaussian quantum 
channels in Sec. In] Then we discuss the solution to the 
capacity for the one-mode case in Sec. |III| and for an 
arbitrary number of modes in Sec. |IV[ Finally, in Sec. [V| 
we analyze the capacity, the optimal states, and the gain 
for Gauss-Markov noise in the full range of the correlation 
strength. The conclusions are provided in Sec. |VI| and 
several mathematical proofs and definitions are provided 
in the Appendices. 



II. CLASSICAL CAPACITY OF QUANTUM 
GAUSSIAN CHANNELS 

In order to transmit classical information via a quan- 
tum channel, one defines an alphabet with letters associ- 
ated with quantum states pj". The quantum channel T is 
a completely positive, trace-preserving linear map acting 
on the input "letter" states: 



(1) 



resulting in output states On average the "letters" 
p'" appear in the transmitted messages with a priori 
probabilities pi so that the overall modulated input state 
is p'" — J2iPiPT- By linearity, the action of T on the 
overall modulated input reads T[p™] = J^iPiPi^*" = 
where p is referred to as the overall modulated output. 
The state p'" is physical only if it has finite energy. 
Therefore it has to obey the energy constraint 



Pi Tr(p-" a^a) < n. 



(2) 



where n is the maximum mean photon number per use 
of the channel, and a and a} are the annihilation and 
creation operators. 

The classical capacity of channel T is the maximal 
amount of classical bits which can be transmitted per 
invocation of the channel via quantum states in the limit 
of an infinite number of channel uses. This quantity can 
be calculated with the help of the so-called one-shot ca- 
pacity, given by the Holevo bound [25] 




P^T[pT]^-Y.P,S{T[pT])^, 



Ci(T) = max 
{pi°>Pi} 

(3) 

with the von Neumann entropy S{p) = — Tr(plogp) 
where log denotes the logarithm to base 2. The max- 
imum in (|3| is taken over all ensembles of {pi^pf^} of 
probability distributions Pi and pure input states p"\ be- 
cause it was proven in [2 6) that the optimal input states 
for noisy quantum channels are pure. 

The term "one-shot" capacity denotes the maximal 
amount of information that can be transmitted by a sin- 
gle use of the channel T. Furthermore, a number of n 



consecutive uses of the channel T can be equivalently 
considered as one use of a parallel n-mode channel T^") . 
Then an upper bound to the capacity of the channel T 
is given by the limit: 



c(r)= lim -Ci(r(")). 

n—^oo n 



(4) 



It has been shown that the latter is the actual capacity 
for particular memoryless O [26l El] and forgetful chan- 
nels [5], but generally is only an upper bound on the 
capacity. Here, we evaluate ([4 ) and find a specific encod- 
ing that realizes this maximal value. For simplicity, we 
refer to C(T) as capacity in the following. 

Let us now consider an n-mode optical channel T^"^. 
In the following, the number of modes of this channel 
corresponds to the number of monomodal channel uses. 
Each mode j is associated with the annihilation and 

respectively, or equivalently 



creation operators aj,a]j. 
to the quadrature operators 



i(aj — aj)l \/2 which obey the canonical commutation re- 
lation \q. 



iSij, where 5ij denotes the Kronecker 5. 



By defining the vector R = (qi, ...,qn;pi, ...,pn) , we can 
express the displacement vector m and covariance matrix 
7 that fully characterize an n-mode Gaussian state p as 



m = Tr [ pR] , 

7 = Tr[(i?-m)p(i?-rn)T]-ij, 



(5) 



J 



/ 

-/ 



where J is the symplectic or commutation matrix with 
the n X n identity matrix /. 

We consider a continuous encoding alphabet, where 
instead of a discrete index we use a complex number. A 
message of length n is encoded in a 2n-dimensional real 
vector 

a ^ (3?{ai},3fi{a2}, ...,3?{a„},9{ai}, ...,3{a„})'''. 
Physically, this corresponds to a displacement of the 
ri-partite Gaussian input state defined by the covariance 
matrix 7i„ (and zero mean) in the phase space by a and 
is denoted by p'^. The Wigner function of p™ reads 



Wt\R) = 



exp [-^{R~ V2a)'^ 7^;^ {R - V2a)] 
(2^)'Vdet (7i„) 



(6) 



where here R G M^" and denotes the coordinates in the 
phase space. In the following we refer to this state as 
quantum input. 

At this point, we follow the standard procedure, 
that is, we only consider Gaussian distributions so 
that the overall modulated input state is a Gaus- 
sian mixture p'" = / d^"'a/(a)p™, where cP'^a = 
d^{ai}d'^{ai} ...d^{an}d'^{an} with (classical) Gaus- 
sian distribution 



/(«) 



exp I 



(27r)"v/det (7,„od)' 
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We refer to the covariance matrix 7niod as classical input 
or classical modulation. We can set the displacement of 
the nonmodulated input state and the mean of the clas- 
sical modulation to without loss of generality, because 
displacements do not change the entropy S{p). Thus, 
since we restrict the set of ensembles over which Ci (T) 
in (|3| is maximized to Gaussian states and distributions, 
we compute the so-called "Gaussian" capacity. This is 
a lower bound to the capacity Ci(T) (hence, for C(T) 
too) which would be the actual capacity if the Gaussian 
minimum output entropy conjecture could be proven to 

hold [21 [in]- 

The action of the channel T*^"-' on an n-mode input 
state carrying the message a reads as in [12, 



out 



rf'"/3/cnv(/3) 



X D{f3n) <E> ... <E> D{/3i) pL" ^^(/^i) ® - ® -D^(/3n), 



(7) 



with /3 = (3?{/3i},...,3fi{/3„},5{/3i},...,3{/3„})T and the 

displacement operator D{(ij) = e^^°'i~^j°-^ . The dis- 
placement is applied according to the (classical) Gaussian 
distribution of the noise /env(/3) with covariance matrix 
7env (which also will be referred to as "environment"). 
If this matrix is not diagonal, then the environment in- 
troduces correlations between the successive uses of the 
channel. These correlations model the memory of the 
channel. 

The covariance matrices of the nonmodulated output 
state and the overall modulated output state p, read, 
respectively, 



Tout Tin ^ Tcnv 
7 = 7out + 7mod- 



(8) 



Note that the covariance matrix does not depend on the 
displacement as seen in definition Therefore, the 

covariance matrix of is the same for all a and the 
output entropy of all displaced states is identical. Thus, 
the one-shot capacity (|3| of this additive channel reduces 
to 



Ci(r("))= sup {5(p) - 5(pr )}• 



(9) 



III. OPTIMIZATION PROBLEM 
A. One mode channel 

In this subsection we consider the case of a single use 
of the channel T^^^ in a similar fashion as in [T^. We 
start our discussion with the results obtained in |16j . 
However, we present the solution for the whole input en- 
ergy domain. We consider the noise covariance matrix to 
have different variances in the quadratures, denoted 7|n^, 
where we choose without loss of generality 7^^^ > 7^^^ 
and off diagonal terms 7^^^ = 0. Any 2x2 covari- 
ance matrix can be reduced to this form by a symplec- 
tic and orthogonal transformation, which changes neither 
the entropy of the output state nor the energy constraint. 
Therefore, by our choice we do not lose generality. As al- 
ready discussed, we restrict ourselves to the optimization 
over Gaussian states. 

In the following we determine the Gaussian capacity 
under the following constraints. The first is the condition 
that p'" is a pure state which together with the definition 
([5]) and the commutation relation imply 



det 7i„ = 



1 



(12) 



The second is the input energy constraint ([2| , that reads 
iL + ll + 7fnod + 7S,od = 2n + 1 ^ A, (13) 

where 7^^^^ , 7^^;^ are the diagonal elements of the matri- 
ces 7in , 7mod and A will be referred to as "input energy" 
in the following. Furthermore, in order for 7?^^^ and 7j^o(j 
to be physical they have to be positive. The optimiza- 
tion problem is solved by using the Lagrange multipliers 
method, with the total Lagrangian being 



- - 3 (^out - \ 
M (7in + 7ta + 7mod + 7mod " ^) ' 



7.9 7P 
' in / in 



(7 



qp\2 



(14) 



with 



We use, in the following, the reduced Holevo ^-quantity 
that reads 



X = S{p) - S{pl 



(10) 



The von Neumann entropy of a Gaussian state p is 
expressed in terms of the symplectic eigenvalues Vj of its 
covariance matrix: 



(11) 



f (a: + l)log(a; + 1) - slogx, , a: > 



We note that for quantum states the symplectic eigen- 
values are always greater or equal to 1/2. 



ill 



qp 



'mod-' ' 



Tout "To u 



(15) 



where 7fjf , 7^^(j denote the off-diagonal terms of matrices 
7in,7modi whcrc 7*'^, 7ou'J denote the diagonal elements 
of 7, 7out and r, p, are Lagrange multipliers. In the fol- 
lowing we summarize the solution of the system of equa- 
tions which correspond to the stationary point of C. The 
details are presented in Appendix We also prove in 
Appendix \K\ that L is concave at this stationary point, 
which implies that the found solution is indeed a local 
maximum of C Though we do not have analytic proof 
that this maximum is global, all our numerical studies 
confirm this. 
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1. Quantum-waterfiUing solution 

First, we find that the solution imphes that the input 
and modulation covariance matrix cross terms vanish, 
that is 7;''jf = 7j^(j = 0. Therefore, the diagonal ele- 
ments 7in^'7mod the eigenvalues of 7i„,7inod- Then 
one obtains the quantum water-filling condition [16j (see 
Appendix A 1 for details) 



7? +7« 
/ 111 ' ' env 



^ ^mod Tin Tt 

which means that the total output energy 

— I /env ' /env 



7mod' 



(16) 



(17) 



is uniformly distributed between the quadratures, that is 



7yH — T^f 



A 
2 



(18) 



We see that this value is also equal to the overall output 
symplectic eigenvalue which will be referred to as water- 
filling level. The optimal quantum input is found to be 
determined by the noise variance ratio, i.e.. 



1 11 



/env 
fcnv 



(19) 



One of the equations links the Lagrange multiplier 



9 i^wf 



= 2^ ^ i7„f = i coth (/i In 2) (20) 



to the water-filling level I'wf, where g' {x) denotes the 
derivative of g{x) with respect to x and In {x) is the nat- 
ural logarithm. In order for the solution to be physical, 
all eigenvalues in ( 16 ) have to be positive. This requires 



an input energy A being above the threshold Athr, i-e.. 



A > Athr — 




(21) 



The optimal eigenvalues which are the solutions of 
(16), (19) are schematically depicted in Fig. [l] (a) for a 
particular noise covariance matrix. The one-shot Gaus- 
sian capacity above the threshold reads 



Ci=g 



/env 



' / en- 



g 



/env Jcnv 



It was shown in [55] that the one-shot Gaussian capacity 
is additive for such a tensor product channel and there- 
fore Ci is the Gaussian capacity C above threshold. 

2. Solution below the threshold 



If the input energy A is below the threshold (21 ), then 
the solution of the Eqs. ( 16 ) and ( 19 ) and the purity con- 



straint ( 12 ) might result in negative modulation eigenval- 
ues, which would be unphysical. Indeed, from the solu- 
tion given by Eqs. ( 16 ) and ( 19 ) one can see that when A 
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FIG. 1. (Color online) Stacked bar plot: optimal eigenvalues 
of both quadratures. Above and below the threshold the op- 
timal input state is squeezed and the less noisy quadrature 
is more modulated, (a) A > Athr. The optimal eigenvalues 
are determined by the quantum water-filling solution. I7„f 
denotes the water-filling level, (b) A < Athr. The modula- 
tion 7^ojj — since 7^^^, > 7?nv ^^wf denotes the "virtual" 
water-filling level. 



decreases the water-filling level I^wf as defined in Eq. ( 18 ) 
[see Fig.[l](a)] also decreases and at A = Athr it crosses the 
level 7? -I- 7|,iv- For lower A, in order to satisfy Eqs. ( 16 ) 



and ([19j) the modulation eigenvalue J^^^ becomes nega- 
tive. In this case the solution of the optimization problem 
lays on the border of the valid physical region. As shown 
in Appendix A 2 we have to set 7^0^^ = [see Fig. [T] (b)] 
and to solve the new optimization problems with a mod- 
ified C taking into account this new condition. Here we 
summarize the results which are presented in details in 
Appendix I A 2| We find that the off-diagonal terms again 



vanish, i.e., 7?f = 7^^^^ = 0. 



Now, one obtains a transcendental equation that solves 
the problem below the threshold, i.e., 



217 



(7^-7^) 



/out q /out 

'ii 



In Appendix |A 3| we derive from this equation that 



(22) 



<7.« < 



/ 7cnv 
/env 



which means that the more noisier quadrature is always 
antisqueezed, and therefore the less noisy quadrature is 
squeezed. We remark that in the limit A — >■ Athr the solu- 
tion of ( p2| coincides with the water-filling solution (16), 
( 19 ). The lower bound is reached for A — > 1 and thus cor- 
responds to vanishing modulation eigenvalues (absence of 
information transmission). We note that Eq. (20) which 
determined the water-filling level reads now 



2V 



(23) 
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This equation now determines fi using the solution of 
( 22 1 . Although we clearly have no longer a water- filling 



solution, we can calculate the quantity z/^f following ( 20 ) 



with /i determined by (23), but refer to it as a "virtual" 
water-filling level. This quantity will have an important 
meaning when evaluating the solution for the multimode 
channel. The reason is that Eqs. ( 20 1 and ( 23 1 in the mul- 



timode channel problem will govern the distribution of 
input energy between the channels, because /i is a mono- 
tonically decreasing function of A. For A > Athr this can 
be seen from the fact that /i is a monotonically decreasing 
function of [see Eq. (20)] and I^^f is a monotonically 
increasing function of A. For A < Athr this is proven in 
Appendix |A 4| This property allows us to relate Athr via 
Eq. ((20l) to 



Mthr = ^g' Q (Athr + 7cV + 7?„v) - ^ j (24) 

such that if A > Athr, then the corresponding /j, < /ithr- 
Moreover, for the lowest input energy A = 1 we can define 
using Eq. ( 23 ) an upper bound /^o for all possible values 




of /i that correspond to A > 1, that reads 



In Appendix |A 3| and |A 4| we draw additional conclu- 
sions from Eqs. ( 22 ) and ( 23 ). We show that 7^ < 7' and 
7^ > 1/2. Furthermore, we find that djf^/dX > which 
means that the antisqueezing is increasing with the in- 
put energy. Moreover, we find that drf/dX > which 
implies together with d'yf^/dX > that the modulation 
in the less noisy quadrature is increasing with A. 

The optimal eigenvalues follow from the solution of 
( 22 ) and are schematically depicted in Fig. [l] (b) for a 
particular chosen noise. The Gaussian capacity below 
the threshold is then calculated by inserting the optimal 
eigenvalues into Eq. (10). 



B. Multimode channel 

We consider in this paper channels with noise correla- 
tions only between the uses of the channel without q — p 
correlations. In the equivalent representation of n suc- 
cessive uses of such a one-mode channel by an n-mode 
parallel channel the covariance matrix of the noise has 
the form 



79 

icnv 





(26) 



where 7^ 



7£ 



are matrices of dimension n x n. The 



absence of q and p correlations is generally considered to 
describe a natural noise. If "f^,^^ and 7^^^ commute, then 
one can diagonalize 7onv via a symplectic and orthogonal 
transformation which changes neither the output entropy 



of the system nor the energy constraint. In the new ba- 
sis the channel becomes a tensor product of monomodal 
Gaussian additive noise channels, for which it was proven 
that the Gaussian capacity is additive [551[3D]. Therefore, 
the optimal input and modulation covariance matrices 
are diagonal, as well as the noise covariance matrix. 

Then the symplectic eigenvalues are functions of the 
eigenvalues of the corresponding covariance matrices: 



■« ^q-.p ^ <!,p , q-.p , q,p 

'III'' II '111,1 ' 'mod,* 'env,z ' 



^q ^p ^q-.p = ^q,p , ^q-p 

'out. 2 'out,z ' 'out.z iin.i ' icnv.i ' 



(27) 



We denote the input energy allocated to channel i as 

= 7in,^ + lfn,^ + ^mod^t + lmod^^■ (28) 

The energy (or photon number) constraint ([2]), (l3^ , for 
n modes can now be written as 



(29) 



The total output energy A of the 71-mode channel is the 
sum of the input energy A and the total energy of the 
noise Aonv (environment) 



X = X + Xe 



i=l 



(7onv,i ~^ 7cnv,i) ' 

(30) 



In jl6j the equipartition of the total output energy be- 
tween the modes was obtained as the solution for such a 
model which we called global water-filling, similar to the 
classical water-filling solution of n parallel classical Gaus- 
sian channels [21]. As in the one-mode case this solution 
only holds above a certain input energy threshold. Now 
we extend the discussion to energies below this threshold. 

The maximum of the Holevo quantity of n parallel 
channels is again determined using the Lagrange mul- 
tipliers method (as in the one-mode case), where we now 
have n purity constraints. 



q p 



(31) 



and the overall input energy constraint, (29). The La- 
grangian is then constructed by a sum of n Holevo x~ 



quantities ( 10 ) for corresponding modes, n Lagrange mul- 
tipliers Ti for the purity constraints, and only one com- 
mon multiplier fi for the input energy constraint. The 
fact that the solution of the system of equations which 
results from the Lagrange multipliers method maximizes 
the Holevo quantity under the given constraints follows 
from the results on convex separable minimization sub- 
ject to bounded variables found in [3T]. This was first 
pointed out for a lossy channel in |17j . In connection 
to our problem it follows that the maximum is attained 
by our solution in the multimode case provided that, for 
the one-mode case, x is a concave function of A on the 
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solution. The proof of the concavity of x in A on the 
one-mode solution is presented in Appendix |A 5[ 

In general, the maximum of the Lagrangian corre- 
sponds to a partition of n modes into three different 
sets, corresponding to one of three types of input energy 
distributions within each mode: the case of a quantum 
water-filling solution with four positive eigenvalues (see 
Appendix III A 1 ), the case of one vanishing modulation 



In this case, Eqs. (16), (19 1, and (34 1 determine the 



global water-filling solution with 



2^ 



and 



eigenvalue with three positive eigenvalues (see Appendix 
III A 2), or the case when both modulation eigenvalues 



q p 

l^out.i — \/ 7onv,i7onv,i 



(37) 



(38) 



vanish and the mode does not participate in information 
transmission (i.e., unmodulated vacuum is sent). We de- 
note the corresponding sets by: A/3, Af2 and Afi. We 
denote the number of modes in the sets as rti, rt2 and 
713, with n = rii + n2 -\- n^. Furthermore, we denote the 
input energies per set by A^^\ A'^\ A^'^-', where 

A(i) = Yl = E ^''^ = E (32) 

ieAfi ieN'2 ii^N's 



If the condition (33) is not satisfied for at least one 



which sum up to the total input energy A ( 29 ) 



mode, then this global water-filling solution has no phys- 
ical meaning because it will lead to negative values of 
some modulation eigenvalues. 



2. Set Ml : Modes excluded from information transmission 

Modes for which both modulation eigenvalues are do 
not contribute to the Holevo quantity or, consequently 
to the information transmission. The zero modulation 



1. Set Ms : Modes with water-filling solution 

For all modes that belong to As the water-filling solu- 
tion described in section UlI A II holds. This means that 
the input energy allocated to each mode cannot be lower 
than its corresponding energy threshold, i.e., 



eigenvalues j^^^^ ^ = 0, i e A/i imply 



^1,? ^ q-.p . 



(39) 



Obviously if the mode is not modulated there is no reason 
to spend input energy on the squeezing of this mode, 
which results in the vacuum state being the optimal input 
state 



A,,; > A, 



thr 



eAA, 



(33) 



where A^ thr reads as in (21) (for all 



Then for all 

I G As the energy equipartition ( 16 ) holds. Moreover 



1 



e AAi. 



(40) 



Eq. ^ guarantees that V^i is the common water-filling ^his is consistent with ^ from which we obtain 



level for all modes due to the common Lagrange multi- 
plier /i, which is a monotonous function of V^i, which 
now reads 



Eq. (40) for A, = 1. 



In order to deduce the set of modes that are excluded 
from information transmission we can use the threshold 



A^ 

2^3^ 



« e A/'s, 



(34) 



r(3) 



where A is the total energy at the output of the modes 
belonging to set A/s. 

As the partition of the input energy between the modes 
is a priori not known the distribution of the modes be- 
tween the sets is also not defined. However, we can deter- 
mine whether a particular mode belongs to set A3 using 
the Lagrange multiplier fi. This is possible, because as 
mentioned before, is a monotonically decreasing func- 
tion of the input energy A and for A^ > Athr.i we have 
fj, < /xthr,i defined in (24) and depends only on the noise 



eigenvalues of mode i. Then we can formalize the defini- 
tion of A3 using /ithr.i as 



A/3 = {j| lithrA > lA- 



(35) 



If Mthr.i > A* for a-U * then the set A3 contains all modes 
and we have 



value /Lto,i defined in Eq. (25 1 which corresponds to Xt = 1 
(vacuum energy) 

A/"! = {i| > A*o,J- 



3. Set M2: Single- quadrature modulated modes 

For the modes for which 1 < A < Athr.j the water-filling 
solution no longer holds. We have to set the modulation 
eigenvalue of the noisier quadrature to in the same way 
as in the one-mode case. Again, as in the one mode 
case we assume that for each mode i the q quadrature is 
noisier than the p quadrature. We can do this without 
loss of generality because, first, for a one-mode channel a 
swap of the noise quadratures does not change the one- 
shot capacity, and second, the one-shot capacity of the 
discussed multimode channel is additive. Then we have 
to set the modulation of the q quadrature for all modes 
belonging to set N2 to 0, i.e.. 



A(3) = A. 



(36) 



7" ^ ■ 

' mod, 2 



0, 



e AAo 



(41) 
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This implies 



(42) 



With the functions /ithr.i and /^o,i defined in (24 1, (25) 
we can simply define this set, i.e., 



(43) 



The eigenvalues that solve the optimization problem 
for the modes of this set are found using Eq. (22), ([23]). 

We note that both, /xo,i and /ithr.z depend only on the 
noise eigenvalues. Therefore, the partition into the three 
sets is completely determined by only one parameter /i. 
Furthermore, we recall that ^ is the common parameter 
which enters the equations for sets M2 , As . 



IV. SOLUTION FOR ARBITRARY NUMBER 
OF MODES 

A. Finite number of modes 

Recall that the solution of the problem for n modes 
is given by the optimal distribution of the input energy 
between the modes. The optimal energy distribution 
within one mode depends on its corresponding set, which 
is given by the noise spectrum and the global parameter 
^. 

Now we present the algorithm that allows us to find the 
solution of our optimization problem. First, we further 
develop the equations that correspond to modes of set 
7V2. We caU the right-hand side of Eq. (122) 



g [i^outA 



2v, 



out, 2 



Tout.z 




(44) 



We note that for given noise eigenvalues / is a function 
of only one independent variable 7fjij. Using definition 



(|28|) and the fact that, for the modes belonging to the 
0, we rewrite 



second set ^'^ 



mod,i 



7?-7f = 



H 



q Q I P 

in.i 7cnv,z ' 7cnv,z' 



Furthermore, from ( |23[ ) we express 
g'iVi - i) ^ 



217, 



71 



(45) 



(46) 



(44) 



Then we insert Eqs. (45) and (46) together into the left 
hand side of Eq. ( 22 ) and obtain with 7^ = 7^^^ ^ and 



Xi between the modes in A/2. Using Eq. (47) and defini 



tion ( 28 ) we can now eliminate variable Xi from Eq. ( 45 ) . 
Thus, we obtain a transcendental equation that deter- 
mines the optimal input eigenvalues jf^ ^ as an implicit 
function of fi: 



7out.iy 1 



= 2/i 



filU 



Now we are ready to calculate the input energies of 
all three sets for a given fj,. First, we evaluate the total 
input energy of "water-filling" modes, i.e. modes in As. 
Using Eqs. (28), (32) and (34) we deduce the total input 



energy used for the modes in A3 as a function of fj. 



7' ■ 

' cnv,2 



7"^ ■) 

' cnv,2 / 



(49) 



Second, the total (vacuum) energy of modes belonging to 
Ai, using (25), reads 



AW(m)= 5] l = ni. 



(50) 



Functions X^^\ii) and A^'^-'(/i) depend on /i through 
A/i,A3 and V^{, which are only functions of and the 
noise eigenvalues. The total input energy for modes in 
A/2 is the sum 



a(^)(m) 



E 



Ai(7in.i'Ai)- 



(51) 



Now we apply the overall input energy constraint 

X'-^^fi) + X^^\fi) + X^^\fi) ^ X. (52) 

Thus, we obtain a closed equation on /i which we can 
solve by iterations. The solution of this system of equa- 
tions provides us with 712 optimal eigenvalues 71^^^ ^ and 
fi which determine all other eigenvalues. Once the opti- 
mal spectra are obtained one can calculate the one-shot 
capacity of n modes {n successive uses) of the channel 
using Eqs. ([9| and (11) 



Ci(r( 



E 

1=1 



g V. 



(53) 



where here i7i,fout,i contain the obtained optimal input 
and modulation spectra. This one-shot capacity of the n- 
mode channel may be considered as n times the capacity 
of the one mode channel. 



HlL,^^ M) = ^^/(7f„.J + 27in,i +7c„v,« - 7on- 



(47) 



Infinite number of modes 



This means that we have established a relation between 
the optimal input eigenvalues 7;^,, i € N2, the global 
parameter /i and the optimal input energy distribution 



In order to make the transition to an infinite number of 
channel uses we have to consider a parallel channel with 
an infinite number of one mode channels, n — > c». In 
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this limit all functions previously labeled with i depend 
now on a continuous parameter x defined on a proper 
domain which depends on the particular model. All sums 
that run from i ~ 1, ...,n now become integrals over the 
whole domain of x. The three sets become now sets of 
continuous variables and cover the whole domain of x; 
they read 



^fl 
^f3 



{x\ no{x) < 

{x\ Hthr{x) < ^ < ^lQix)}, 
{x\ fithr{x) > fJ,}, 



(54) 



where ^thiix)^ fJ.o{x) are defined as in (24), (25 1 where 



index i is replaced by x. Equations (48)-(51) remain the 
same, except for the replacements 7in j by (x) and the 
sums over i by integrals over x. 

Once the /i is found which is the solution of ( |52[ ) we can 
determine the optimal spectra 7in^(a;) and 7i'^^dT^)- The 
found optimal spectra are used to evaluate the capacity 

C= lim -Ci(r(")) 

_ 1 

(55) 

where A is the spectral domain of x and \A\ is its size. 



dx { Q { i^{x) - 1:] - 9 [ i^out{x) - 



In the case of a global water-filling, i.e 
then ( 55 ) simplifies with ( 37 1 and dSS" 



if fithri^) > Va; 
to 



C^g 



1 



2\A\ 



-7^ 
/ cnv 



(^)} 



xeA 



1 



dxg y\J Jcnv{x)jinvix) 

where "tcrivix) ai'e the noise eigenvalue spectra. 



(56) 



V. GAUSS-MARKOV CHANNEL 

Now we consider a particular channel where the cor- 
relations of the noise are modeled by a Gauss-Markov 
process. Note, that for the whole class of noises with cor- 
relations given by stationary (shift-invariant) Gauss pro- 
cesses, the covariance matrix that describes such noise 
is a symmetric Toeplitz matrix. A straightforward ex- 
ample is the Markov process of order V, also called au- 
toregressive process with white Gaussian noise [55] (see 
Appendix [C]). Our above treatment applies to such chan- 
nels. However, in the following, we restrict to the sim- 
plest case of ~ 1. We have already introduced the 
covariance matrix of such a noise in ,16J and state here 
only its definition. Essentially this is a classical noise 
with nearest-neighbor correlations in the q quadratures 
(anti-correlations in the p quadratures). Its covariance 
matrix contains Toeplitz matrices for the q quadrature 
and p quadrature, so that 



_/Af(0) 
^''"^ ~ \ M{-4>) J ' 



(57) 




(b) 

e 









7env(*^ 













IT/2 

X 



FIG. 2. (Color online) Stacked area plot: optimal input and 
modulation eigenvalue spectra y^^{x), "f^^di^) noise spec- 
trum 7lnv(a;) (for p in dashed and dotted curves) for a par- 
ticular (j) and N. (a) Global quantum water-filling solution, 
where A > Athr. J^wf (solid bar) denotes the water-filling level, 
(b) Below threshold: Modes with x £ [a, tt — a] belong to As 
with water-filling level I7„f (solid bar). Modes with x £ [0, a] 
and X £ [n — a,n] belong to set M2, where the modulation is 
below i7„f. 




FIG. 3. Functions ij.o{x) (upper curve) //thr(a;) (lower curve) 
for (a) (j> = 0.5, (b) = 0.7, (c) (j> = 0.9, (d) (j) = 0.99. For all 
plots we took A = 1. 



with 



0<d)<l, N>0, 



(58) 



where is the correlation parameter and N is the vari- 
ance of the noise. We remark that a Toeplitz matrix has 
the same values on each k = \i ~ j\th diagonal. 



Matrix (57) can be diagonalized in the limit of 
n ^ 00 using a passive symplectic transformation |16j . 
which allows us to study the channel system in the diag- 
onal, noncorrelated basis, since entropies and the input 
energy constraint remain unchanged by such transforma- 
tions. 
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The spectra of the noise quadratures in the hmit of an 



infinite number of uses of (57 1 read 



1 + 2(l)cos{x) ' 



e[0,2^], (59) 



with the upper sign for the q quadrature and the lower 
sign for the p quadrature. As this spectrum is mirror 
symmetric with respect to x = tt and since the Gaussian 
capacity of the channel is additive, we reduce the spectral 
domain to a; G [0, tt]. 

In order to find the capacity of the channel for given 
noise parameters (j), we first check whether the thresh- 
old condition ( |33[ ) (here with continuous parameter x re- 
placing index is fulfilled for all x. If this is the case, 
then the solution is a global water-filling, depicted in 
Fig.[2](a), where optimal eigenvalue spectra are obtained 
by ( |l9| , (16) (for all a;), and the capacity (56) can be 
simplified to [IS] 



C = g(n + N)--J dxgy\J jcnv {x)'^cnv (x) 



n > 



2<j) 



N- 



(60) 



1.5 



e 



0.5 - 



.-'2- 



uJ2 
X 



FIG. 4. Functions /ithr(a;) (lower dashed curve), fio{x) (upper 
dashed curve) and values for /i (solid bars) for different input 
energies A, and noise parameters = 0.85, A'^ — 1. From 
top to bottom the values are p — 1.45(A = 1.006), 1.34(A = 
1.04),0.42(A = 3),0.04(A = 35). The numbers indicate the 
intervals on the a;-axis that belong to sets Ni,N2 or As. 



If the threshold condition is violated then we can apply 
the algorithm which is presented in section [TV] 

We note that for different noise parameters the thresh- 
old functions Ho{x), ^ith^ix) may have a complicated pro- 
file as depicted in Fig. |3] for different noise parameter 
values. In Fig. |4] we illustrate, for a particular choice 
of the noise parameters (iV, 0) different partitions of the 
spectral domain between the sets for different input en- 
ergies A corresponding to different fi. 

Our result confirms that the modes belonging to Af2 
are squeezed in the less noisy quadrature which is the 
one that is modulated, as depicted in Fig. [5] An ex- 
ample plot of optimal input and modulation spectra for 
A < Athr is shown in Fig. [2] (b). Wc see the naturally 
expected behavior of the capacity in Fig. |6] It decreases 
with increasing noise variance iV and increases with in- 
creasing noise correlations 4>. We note that the capac- 
ity increases with (p up to the noiseless capacity at "full 
correlations" (</)—)■ 1). This limit will be discussed in 
section IV Bl 



(a) 



nil 
X 




FIG. 5. Optimal input 7in(a;) (solid curve) and modulation 
7mod(^) (dashed curve) eigenvalue spectra of the q quadrature 
(p quadrature spectra are the same but mirrored with respect 
to a vertical line at 7r/2) vs. spectral parameter x, for = 
0.85, A/" = 1 and A < Athr. The partitioning in sets is taken 
from Fig. |4] In (a): A = 3 which corresponds to /i = 0.42, in 
(b): A = 1.04 which corresponds to /x = 1.34. 



fcth diagonal reads 



A. Optimal quantum input state 



'in,fc 



27r 



t^a^e* 7f/(a;), fc = 0, 1, 2, oo. (61) 



An important issue is to derive the covariance matrix 
of the optimal input state in the original "correlated" ba- 
sis. We know that in the basis where the noise, modula- 
tion, and input matrices are diagonal, the optimal input 
spectrum corresponds to a product of one-mode squeezed 
states. By using general properties of Toepltiz matrices 
we conclude in Appendix[B]that in the limit of an infinite 
number of modes the optimal input covariance matrix in 
the original basis is also Toeplitz and we found that its 



We can express ^'{^{x) exactly in the case of global water- 
filling (A > Athr) which we consider for the rest of this 
subsection, that is, 



hcnvix) 

7™v (x) 



Inserting the definition of the noise spectrum of the 



Gauss-Markov channel (59) we deduce the spectra for 
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FIG. 6. (a) Capacity C (in bits) vs. correlation where 
from top to bottom A'^ = 1, 2 and 3. (b) Capacity C (in bits) 
vs. noise variance A'', where from top to bottom (j) = 0.9,0.7 
and 0.5. The input energy is A = 3 for both plots. The 
dashed part of the curves corresponds to the global water- 
filling solution with A > Athr- One observes that the capacity 
for full correlations — >■ 1 tends to the capacity of the ideal 
noiseless channel TV = 0. 



the q and p quadrature of the optimal input matrix, i.e.. 




(62) 



where the upper sign is for q and the lower for p. In order 
to verify that the overall state is entangled we can check 
whether the reduced single mode states are mixed, i.e. 
whether for the reduced covariance matrix we have 



det 7i, 



q p 

~ 7in,0'^in,0 



> 



1 



> 0. 



(63) 



Integration over the whole domain to 27r leads to 
iLfi = 7f„,o; Then we find for 7fn,o(<^ = 0) = 1/2, which 
means that in the absence of correlations the optimal in- 
put state is a set of coherent states and not entangled. 
The limit of </> — ^ 1 is unphysical because in this limit 
each single mode state becomes a thermal state with its 
temperature tending to infinity. This corresponds to an 
overall maximally entangled state. It is easy to show that 
( 63 1 is monotonically increasing from <j) — Q to cf) = 1 and 
therefore we conclude that for all > the optimal input 
state is entangled. 

In order to express the covariance matrix of the overall 
modulated output, let us recall that in the global water- 
filling case the overall modulated output eigenvalues are 
identical [7''^(a;) is constant in x\ and we can express 
them using Eq. ([37]) as ^'^'^{x) = n + N + 1/2. Therefore, 



from Eq. (61) we easily see that only the main diago- 
nal (A; = 0) has nonvanishing values and these values 
are identical. Then the covariance matrix 7 is propor- 
tional to the identity matrix /, and therefore it is di- 
agonal in the initial as well as in the rotated basis is 
^I'P ^ I{n + N + 1/2). This means that the sum of the 
optimal input and modulation covariance matrix has to 
cancel the correlations of the noise. 



B. Full correlations 

We observe in Fig. [6] that for fixed N and A the higher 
the correlations are, the higher is the capacity. Further- 
more, for — >■ 1 the capacity tends to the capacity of the 
noiseless channel 



lim C = Co 

4,^1 



9{n), 



(64) 



where the noiseless capacity Co was calculated in e.g. 
Equation ( 64 ) can be deduced by the following reasoning. 
For any < (f> < 1 the capacity C is upper bounded 
by Cq. In addition, C is lower bounded by the optimal 
transmission rate when using coherent states, which is 
denoted by R in the following. Thus, we need to show 
that for — > 1, both bounds fall together. 

R is easily calculated, as the restriction to coherent 
input states basically leads to the discussion of two clas- 
sical multimode Gaussian channels with noises "f^^{x) = 
7cnv(^) + 1/2- Clearly, the solution of the optimiza- 
tion problem completely reduces now to the classical 
water-filling [24J which determines the optimal modula- 
tion spectrum. For the noise spectrum of the Gauss- 
Markov channel (59) we find the global water-filling so- 
lution 



7r 

R = g{n + N) - ^ J dx g (^^J jenv{x)jinv{x) - ^ 



26 

n> -—N. 



(65) 



Since for global water-filling the overall modulated out- 
put state is identical in R and C ( 60 ) the difference to the 



capacity comes from the difference in the nonmodulated 
output only and is remarkably little. Indeed if we look 



at the output eigenvalue spectrum for R in ( 65 1 which 
reads 



i'outCa;) = 
instead of 



/l 7cnv(a;) +7cnv(a;) q / \ p / \ 

y 1 + 2 7onv(a;)7cnv(2:) 



for C ( 60 ) we see that the two formulas simply differ by 



the terms which are the arithmetic mean of the noise 
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eigenvalues 7|„^(a;), 'ji^^ix) in the first one and the geo- 
metric mean in the second one. As the geometric mean 
is always less or equal than the arithmetic mean one con- 
firms that C > R. 

Below the energy threshold one has to solve 



-7env(a) 



1 



(66) 



for a (depicted in Fig.[2](b) but here with ^f^^{x) = 1/2) 
which is the a: value that defines the sets Af2 and As when 
restricted to coherent states. For the found value of a we 
determine the optimal modulation eigenvalues 



7mod(2;) = 

7mod(a;) 



u[X - u;[7p-„^(_ct; - ' /onv V-^^J i 

9{n-a- x)[7«„^(a) - 7o„v(a;)], 



")[7cnv(a) - 



7onv i"^)] 1 



(67) 



where 0{x) is the Heaviside step function. By inserting 
lini^) = l/2,7mod(^) and ^l£,{x) in ^ one obtains R 
for n< 2N(j)/{l- cj)). 

In the limit of full correlations — > 1 the noise spectra 
7|^^(a;) tend to for < x < vr and to infinity for x = 
(for the g-spectrum) and x = tt (for the p-spectrum) . Due 
to the finite energy of the noise, ^ dxj^mi^) = 
these functions become delta-like distributions. In this 
limit a — >■ and the solution to R is given by a classical 
water-filling over a vacuum noise spectrum with infinite 
edges which though can be shown to give an infinitesi- 
mally small contribution and therefore can be neglected. 
Thus (64) is proven. The same result was obtained in 



[33] for a channel with additive Markov noise, which be- 
comes a collection of thermal channels when the noise is 
diagonalized. 



C. How useful are the optimal input states? 

In this subsection we evaluate the gain G from the use 
of the optimal input states with respect to coherent prod- 
uct states for the Gauss-Markov channel for two modes 
and an infinite number of modes. Our motivation here is 
that the optimal input states are entangled and therefore 
may be not easy to generate. On the contrary, coherent 
states are easily accessible in the laboratory by standard 
tools of quantum optics. The gain G is given by the ratio 
of the capacity C to the optimal transmission rate using 
coherent states R (discussed in section V B I 



G = 



C 
R' 



(68) 



The gain was already discussed in 12J for the case of 
a two mode additive channel which is identical to our 



Gauss-Markov channel with noise covariance matrix (57) 



taking n = 2. We remark here, that the capacity of two 
modes with correlations (f> (in q and p, but no q — p cor- 
relations) is identical to the capacity of the mono-modal 
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n 



FIG. 7. Gain G vs. n for an infinite number of modes and 
for two modes (inset). For both plots we took SNR = 3 and 
(f) — 0.7 (solid curve), 4> — 0.9 (dashed curve), (f> = 0.99 
(dotted curve). 



phase dependent channel discussed in section |111 A[ be- 
cause in the two-mode case the diagonalized noise spec- 
trum leads to two phase-dependent monomodal channels 
with inverse variance in q and p. As it was shown '12^ the 
gain for such a channel exhibits a maximum with respect 
to the input energy constraint n for fixed signal-to-noise 
ratio 



SNR = 



n 

N 



and correlation (p. Furthermore, it was deduced that the 
gain increases with increasing correlations (/> between the 
two modes. 

In the case of an infinite number of modes, we know 
already that in the absence of correlations the optimal in- 
put states are coherent states, and therefore there is no 
gain {G — 1). For full correlations, the behavior is essen- 
tially different from the two-mode case: since the chan- 
nel becomes effectively noiseless, coherent input states 
are optimal in this limit as well, whereas for two modes 
the highest available squeezing is best. Therefore, an in- 
teresting question is where we find the maximum gain 
with respect to the noise correlations in the limit of an 
infinite number of uses of the channel. In Fig. [T] we 
plotted the gain G vs. n for fixed SNR and different 
for an infinite number of modes and for two modes. We 
see that unlike in the case of two modes, where the gain 
with higher correlations is always higher, in the case of 
an infinite number of modes the maximum of the gain 
is found for some intermediate correlations. However, in 
this plot one does not see the dependency on the SNR. So 
the question that follows is: What is the dependence of 
the maximal gain (with respect to n) on (f) and the SNR? 
In order to answer this question we make a contour plot 
of maxnG vs. (p and SNR. In the case of two modes we 
see in Fig.|8]that the optimal gain is obtained at ^ = 1 for 
a certain SNR. In addition, in this case, the increase in 
gain at high correlations is very strong compared to that 
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FIG. 8. (Color online) Contour plot of the maximal gain 
max^iG vs. 0, SNR for two modes. In the area above the 
dotted line the quantum water-filling solution holds. 
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FIG. 9. (Color online) Contour plot of the maximal gain 
max^i G vs. <\>, SNR for an infinite number of modes. In the 
area left to the dotted line the global quantum water-filling 
solution holds. 



VI. CONCLUSIONS 

We have presented an algorithm for calculation of the 
classical Gaussian capacity of the Gaussian channel with 
additive correlated noise. This method is applicable to 
all cases where there are no q — p correlations in the 
noise covariance matrix and, moreover, the q and p blocks 
commute at least asymptotically in the limit of an infinite 
number of uses. This applies, in particular, to the whole 
class of channels in which noise correlations are given by 
a stationary Gauss process. 

We applied this method to a channel with a Gauss- 
Markov noise as defined in [TB] , which has asymptotically 
commuting block matrices. We found that in the limit of 
full correlations the capacity tends to the noiseless capac- 
ity. We calculated the covariance matrix of the optimal 
input state not only in the eigenbasis of the noise co- 
variance matrix but also in the original, correlated basis. 
In this correlated basis the optimal input state is entan- 
gled and we found that the degree of entanglement scales 
with the correlation parameter of the noise from no en- 
tanglement (i.e., a set of coherent states) to a maximally 
entangled state. 

Furthermore, we discussed the gain from using optimal 
entangled input states with respect to coherent product 
states in the case of two modes and an infinite number of 
modes. We found that, contrarily to the two-mode case, 
where the gain always strongly increases with correlations 
for any SNR, for an infinite number of modes a high gain 
is achieved in a region of high correlations and low SNR 
and of lower correlations and higher SNR. In addition, 
the gain in the limit of an infinite number of modes does 
not drop as sharply with decreasing correlations as in 
the two-mode case. We also observed that a Gaussian 
modulated coherent-state encoding already achieves not 
less than 90% of the Gaussian capacity. 



at lower correlations. Furthermore, a low SNR seems to 
benefit from entanglement more than a higher SNR. 

For an infinite number of modes the situation is dif- 
ferent, as we can see in Fig. [9} instead of a sharp edge 
toward high correlations we see an almost-flat area of 
maximal gain in the region of high correlations and low 
SNR. This holds, on one hand, for a high correlation and 
low SNR but, on the other hand, also for less correlated 
noise and a higher SNR. Furthermore, the enhancement 
is rather robust and does not drop as sharply with de- 
creasing correlations as in the two- mode case. However, 
as the region of high gain has input energies below the 
global water-filling threshold, where the optimal input 
squeezing becomes quite complex [as depicted, e.g., in 
Fig.[5](b)] a modulation of coherent states might be prac- 
tically more favorable, because it is already quite efficient 
as the gain due to entanglement does not exceed 10%. 
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Appendix A: Results for one mode 

In the following we present the solution via the La- 
grange multipliers method for the optimization problem 
for the one-mode channel introduced in section ITlll 
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Search for the extremum 



The extremum of the Lagrangian £ defined in ( 14 ) 
must satisfy 



V£ = 0, 



where 



f d d 



d 



d 



d 



d 



This corresponds to six equations: 



qp 

mod 



K.{v) -f - K(l^out)7out 



/in 

v9 - 



T7r„ = 0, (Al) 



<V) 7« - «(i.out) 7o'ut - M - T7f„ = 0, (A2) 

K(l7)7P-/i = 0, (A3) 

k(I7)7''-M = 0, (A4) 

-'^l^) (7^ + 7rod) + '^(^out) iZ + = 0, (A5) 

-«(^)(7r + 7rod) = 0, (A6) 

where r are Lagrange multiphers and 



2x 



From Eqs. (A3) and (A4) we derive Eqs. (16 1 and (20 1 



Since k{x) > for all x > 1/2 we find from Eq. (A6) that 
-fZ — ~7nfod- Therefore, Eq. (A5) simplifies to 



7i'5f[«:Kut) + r]=0. 



(A7) 



If one assumes that 'jZ then the resulting multiplier 
T leads to a contradiction when inserted in Eqs. ( Al ) and 



(A2). Thus, we conclude that 7*f = 7;^^^^ = 0. Finally, 



by combining equations (A1)-(A3) one deduces Eq. (19 1 



In order for the solutions to be physical the modula- 
tion eigenvalues 7^od''>'mod have to be positive, which is 
the case for an input energy above the energy threshold 
(21). For such A we will show now that £ is concave 



on the solution which proves that we found indeed a lo- 
cal maximum. For fixed A, we can prove concavity by 
checking whether the Hessian H{C), i.e. the 6x6 matrix 
that contains the second derivates of C with respect to 
all variables 7?^;^, 7^^^^, 7in^mod' negative definite. First, 
we find that 



0, 



in, mod, 



which correspond to the derivatives of the left hand sides 



of Eqs. (A5) and (A6) with respect to 7i„ ,7^od- There- 



fore, we cab write the Hessian H{C) in the block form 



H{C) 



i/cov(/:) 



where H^g_^{C) contains only second derivatives with re- 
spect to 7in^,7mod ^"^"^ Hcov{C) second derivatives with 



respect to 7, 



in. mod ' 



By using 7,"^^ ^ 



Eqs. (Al), (A2| we find 



A + B A 
A A 



where A = K{iy), B — K(j/out)c, with 

' env ' / e 



L- — '^y J env Ji 



env q p 
^env '7env 



The eigenvalues of ( A8 1 read 



hi 



IP _ 



-A--± \IA^ + 



52 



and 



(A8) 



(A9) 



(AlO) 



which are both negative, since A,B>0. 

Now we show that the eigenvalues of TJvar {^) are also 
negative. However, instead of considering the Hessian 
i/var('C) we consider equivalently the Hessian _ffvai(x), 
where we embedded constraints (12), (13) in which 
becomes then a function of only two variables 7inj7^od 
and thus iJvar(x) is a 2 x 2 matrix. Then we find that 
-ffvar(x) has the same shape as (A8), where now 



A^ 



9' {V 



B = 



g' (i^out 



/env 



^out4(7fJ3 



(All) 



Since A,B > Q we conclude again that both eigenvalues 
[that read like (AlO)] are negative. Therefore, the total 
Hessian H{C) on the solution is negative definite, which 
proves the concavity of C at the extremal point for an 
input energy above the threshold ( 21 ) . Thus, we conclude 



that we have found a local maximum of C. 



Below the threshold 



If the input energy is below the threshold (21 ) then the 



extremum of the Lagrangian lays outside the physical re- 
gion. In this case the maximum lays on its boundary 
which corresponds to one or both modulation eigenval- 
ues being 0. If both modulation eigenvalues are no 
information is transmitted which is clearly not the opti- 
mum for an input energy A > 1. Now we put 7.' ^ — 



which replaces Eq. ( A3 ) (if we put instead 7^^^ = we 



would obtain in the following the same set of equations 
up to a permutation of indexes p and q). Therefore, this 
degree of freedom does no longer exist and there is no 
q — p correlation in the modulation: 7^^ = 0. Now, 
we have to find the extremum of the Lagrangian again, 
where the new gradient reads 



V = 



d d 



d 



d 



dlL ' ^7^ ' 57mod ' ^7* 



The previous equation, (A5), is now simplified to 
-lZ\^i^) - >^{vout) - r] = 0. 
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Again, if one takes 7^ ^ then the solution for r when 
inserted into Eqs. (|A1[) and ( A2) leads to a contradiction, 



namely 7^^^^^ < 0, which is unphysical. Therefore, we 
conclude that 



7.''^ = 0. 

/ m 



(A12) 



and h decreases. Therefore, by doing this we increased x 
which leads to a contradiction because we the assumed 
optimal solution is in fact not optimal. Thus, for 7'' > 7^ 
the lemma is proven. 

Now we consider 7' < 7^ and assume that wc have 
found an optimal solution. Then from Eq. (A13) it fol- 
lows that 



By inserting (A4| into ( Al ) one obtains t which when 



inserted into ( A2 1 leads to the transcendental equation 



( 22 1 . The Lagrange multiplier [i is given by ( A4 ) leading 



to Eq. (23) 



7- > 

/in ^ 



1 



'To 



(A16) 



From Eq. ( 22 ) we can deduce a bound on the optimal and from the condition 

ily that 7f„ > 1/2 > 



input squeezing. From the fact that 



/cnv 



> 7^j-,y we deduce eas- 
Taking into account that 
7mod > 7mod = WC coucludc that in fact 7*^ > 7^ which 
contradicts to our assumption. 

Thus, the lemma is proven. □ 



< 



as 1/x g' {x — 1/2) is a monotonically decreasing function 
and V > i^out, it follows that 



A direct consequence of Lemma [T] and Eq. (A13) is 
that 



|7^-7^|>l7Lv 



7' 

/cnv 



(A13) 



where the equal sign holds if both left and right hand 
side are and the expressions inside the absolute values 
on both sides have the same sign. When A < Athr t hey 
cannot be 0, because otherwise 7?^ = l/2y/^env/jinv and 
TyP — Tyi leading to A = Athr which contradicts to our 
assumption that we are below the threshold. This proves 
that 



7Vr 



(AM) 



Let us remind that up to here the obtained results are 
invariant (up to the permutation of indices q and p) with 
respect to the choice which modulation eigenvalue is set 
to 0. 

Now we prove that the choice which modulation eigen- 
value has to be set to is determined by the relation 
between jg^^ and 7^^. 



Lemma 1. For an input energy A < Athr under the con- 
dition 7|jjy > ^P^^ the maximum of x given by Eq. ( 10 ) 
is achieved when 



7Ld = 0- 



(A15) 



Proof. Let us assume that in contrary to the statement 
of the lemma we have ^^^.d = 0- We know that due 




in(thr) 



(A17) 



meaning that the squeezing of the input state for A < Athr 
is always smaller than for A > Athr- 

The next step is to prove that the extremum we found 
is a maximum. In order to do this we verify that C is con- 



cave for an input energy below the threshold (21 1, where 
A is fixed. First, we see again that all cross derivatives 
with respect to and the three variances vanish on 
the solution. Therefore, we consider in the Hessian H{£) 
the second derivative with respect to 7?^ separately. By 
using Eqs. (A2| and (A4) we find 



9(7^)2 



1 < 0, (A18) 



since 7lt > 7in- 

Again, for the remaining three variables instead of con- 
sidering the Hessian of (£) we consider equivalently the 
Hessian of x with constraints ( 12 ) and ( 13 ) embedded. 
Now, X becomes a function of only one variable 7?^, and 
thus the Hessian is only the second derivative of x with 
respect to j?^. We prove in Appendix 



A4 



that 



d'x 



dF 



< 0, 



(A19) 



to Eq. (|Al4p we have two possible cases: 7"^ > 7^ or where F is given by (|A3lj). Therefore, the full Hessian of 



j<i < jP, 

Let us assume first that 7'' > 7^ and that we have 
found an optimal solution. Suppose we remove a frac- 
tion of 7^o(j which is smaller then half of the difference 
7'' — 7^ and set j^^^^ equal to this fraction. This will not 
change the input energy as well as the output entropy 
(second term in x). However, the overall modulated out- 
put entropy (first term in x) will increase. The reason 
is that for a constant arithmetic mean of a and b the 
geometric mean increases when the difference between a 



C is negative definite for input energies below the thresh- 
old as well. Thus, we have shown that C is concave on 
the solution, which proves that we indeed found a local 
maximum of C. 



3. Bounds below the threshold 

In the following we deduce bounds on the overall mod- 
ulated output variances for A < Athr- From Eqs. (A13) 
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and (A17l it follows directly that 
7^ <7«. 



(A20) 



Furthermore, we can find a lower bound on 7^. Sup- 
pose < 1/2, then we have 



1 
V 



1 



9 V 



> 



>9' 



1 




2 ' 2 



(A21) 



since g'{x — 1/2) is a monotonically decreasing function 
of X. Thus, we can rewrite Eq. (22), 



9 ['^ 



U ll^q \ 2 

2 ' 



(A22) 

and by using (|22| and the fact that below the threshold 



f = 7out we find 



2l'out 

where 



E > 



2 /out 



Tout 



^9 

/cnv 



7c" 



(A23) 



(A24) 



Thus, our assumption 7^^ < 1/2 leads to an inequal- 
ity which depends solel y on 7^^ , 7^^^ , 7^^^ with the con- 
straints on 7!^^^ given by (A17) and for the noise variances 



that 7?nv > 7?nv, < 7?„v < 1/2 - l/(47i'^J. Since, 
inequality ( A23 ) is always violated for the given con- 



straints, we come to a contradiction which proves that 
7P > 1/2. (A25) 

4. Monotonicity of fj, 

Lemma 2. The Lagrange multiplier ^ is a monotonically 
decreasing function 0/ the input energy A on the solution, 
and moreover 



dfi 
dX 



< 0. 



(A26) 



Proof. For A > Athr the proof follows directly from the 



definition of 11 by Eq. ( 20 ) , the definition of v^jf by 



Eq. (18 1 and the fact that g'{x) is a monotonically de- 
creasing function. 



Now we prove the lemma For A < Athr- Using Eq. (23) 
we can write 



We can upper bound this quantity if we use the following 
property of 5(2;): 



g' (17 -1/2) 



>.g"(l7-l/2). 



This leads to 



dx 2173 ' dX' 



(A28) 



(A29) 



Since all factors in (A29) except for d^^/dX are clearly 
positive the lemma will be proven if we show that 



(A30) 

This derivative can be expressed in terms of the function 



217 



2t'out 



p 

out g /out 



Indeed, when Eq. ( 22 ) holds then we have 



dX 



OF 
dX 
dF 
07' 



(A31) 
(A32) 



We observe that 



dX 4172 



7^7^-7') 



(A33) 



because g"{x) < 0, g'{x) > and j'' > 7^'. Thus, in 
order to prove inequality Eq. (A30) it suffices to prove 
that 



dF OF 



dX 97?„ 



< 0. 



(A34) 



By carrying out the partial derivatives we rewrite 
Eq. (|A34|) in the form 



-^7^(7-7^m-^T.<0, 
where Ti and T2 denote the expressions 

1 



(A35) 



(A36) 



+ g Vont 



(7i 



q \3 ''out 



and 
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Observe that all factors in front of Ti and T2 are positive, 
since 7^, J?, t'out,7' — 7^ > 0- If we prove positive of Ti 
and T2, then inequality (A35) will be proven as well as 
the lemma. 

The positivity of Ti can be verified via its partial 
derivatives with respect to 7^,7'' which lead to: 



dTi 



= {7" +7"[3 - 4r(r + 37'')]}Tn(7^7^), (A38) 



(7'-7^)Ti2(7^7^), 



(A39) 



where 



rii(7',7'') 
^12(7', 7^) 



(7«+ 7^)3(1 -47^)^ 
4(7'')^(l + 4(7'-)^) 

VfV(7' + 7^)^(1 - 47^)^ ' 



(A40) 



Clearly the two functions Tnj^'^ ,^ p) an d ri2(7'?,7^') are 
positive for all 7'''^ > 0. Then (A39) is also positive 



since 7' — 7^ > 0. From equations (A20) and (A25) we 



have that > 1/2, 7« > 1/2. Then, for all these values 
of 7', 7^ it is easy to v erify that the factor in front of 
I'ii(7'',7^) in Eq. ([AM]) is negative. Thus, 



dTi/dj'^ < 0, dTi/d^P > 



(A41) 



and Ti takes its minimal value at the boundary of the 
allowed region for 7^,7^, namely, at the point where 7'^ 
is maximal and 7^ is minimal. Observe that for M2 using 
Eqs. (|A17|, (|A20P and (|A25| we have 



Now we show that T2 > as well. We rewrite T2 as 

,2 



=C 



{z' - l)^/3(z^out) + 



(A44) 



where 



7in(thr) 

g" {fput - ^) 
.9' (i^out - I) 



9' (^out - (7eV)' 



(A45) 



where 7f,j(^jjj.') was defined in Eq. (A17) and therefore 



< z < 1. Using these notations we express 
Vout = Vl/4 - i^cnv/2{z + 1/2) + 



(A46) 



and since ^ > we can rewrite the desired inequality 
r2 > in the equivalent form 



Kz, Vout) > 



where 



/3(t'out) 



32z2 



(A47) 



- Z2)2 + 16z2l/2^j - 1 - z2)(l _ ^2)2 

(A48) 

We observe that 



lim h{z, Vo^t) 



16 



^-^0 ~ 4^^,t-l 
It is easy to check that the inequality 



16 



> 



l3{Vont) 



(A49) 



(A50) 



1 



' /cnv 
/cnv 



+ 7, 



9 

cnv • 



(A42) 



Then Ti takes its minimal value for 7^jj„ — 1/2 and 
7?nax = 5\/7cnv/7?nv+7lnv Therefore, if Ti at this point 
is positive for all values of 7|„v7 7cnv then it is positive in 
the whole allowed region of 7* and 7^ . 

In order to evaluate Ti at this point we derive it with 



respect to 7, 



env ' 7cnv 



97, 



and we find that 



cnv 



(A43) 



u /env 



holds Vi^out > 1/2. Thus, if dh{z, Vant) I dz > holds 
for all fout and z in the allowed region then the desired 
inequality (|A47|) holds. We find that 



— 64z- 



dz 



b{z, I'out) ' 



(A51) 



where 



2z'' 



a{z, i^out) = -1 - 82^1/2^^ - 32''(8z^^ut - 
+ ?(z,i^out)(l + z2 + 2z4), 

&(^, J^out) = [z^ - I^out)(l + Z^ - l{z, Vout))'^, 



1{Z, Vont) = \Jl + Z^ + 



2z2(8,y2 



(A52) 



Then, again Ti takes its minimal value at the point where 
^cnv is maximal and ^^^^ is minimal. At this limit point 
we find 



>00,7c: 



0, 



where the limit is reached from above. This proves that 
Ti > 0. 



Clearly I'out) is negative and therefore, if a(z, i-'out) is 
negative as well then dh {z, Vo ut)/dz > 0. Since the first 
line in 0(2, fout) in Eq- (A52) is negative in the allowed 
region of ^^out and 2, and the second line is positive we can 
make a comparison of squares of the first and second line, 
which confirms that indeed 0(2, Vput ) < 0. Thus, T2 > 0, 
as well as Ti > which proves (A34) which means that 
(A30I holds and thus, the Lemma is proven. □ 
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Additionally by combining Eqs. 
( A34 ) we conclude that 



(A32), (A33I and 



dX 



> 0. 



(A53) 



This means that the antisqueezing in the more noisy 
quadrature is always increasing until the squeezing value 
at A = Athr is reached. 



5. Concavity of the Holevo x-quantity in A 



Lemma 3. The Holevo x-Quantity given by Eq. (10) is a 



concave function of \, on the solution of the optimization 
problem. 



Proof. For A > Athr we find that 'j?^ is given by ( 19 ) and 



independent of A. Therefore, we conclude from Eq. ( 18 ) 



that X is on the solution a function of only one variable 
A. Then, at the extremum of C the second partial deriva- 
tive of X with respect to A is equal to the total second 
derivative, which reads 



9A2 



< 0. 



(A54) 



Thus, we have shown that above the threshold x is a- 
concave function of A. 

For an input energy A below the threshold, depends 
on A via the implicit function given by Eq. (22). There- 



fore, the total second derivative of x with respect to A 
has to take into account this dependence. Now this reads 



dA2 



(Px 
9A2 



rf7in , dx 



9A97?„ 
d^X . 



92 X 



" 'in 



dX^ 
dX I dX 



(A55) 



One can easily show using Eq. (A4) that dx/dX — fi 
and by E q. ^KM\ it follows that 9x/a7f„ = F. Thus, 



Eq. (A55) simplifies on the solution to 



d^ 
dX^ 



dii 
dX 



< 0, 



(A56) 



as proven in Appendix | A 4[ which proves the lemma. □ 



Appendix B: Eigenvectors and Eigenvalues of 
Toeplitz matrices 

In the following we derive the optimal input covariance 
matrix in the case of global water-filling. 

All Toeplitz matrices that belong to the Wiener class 
commute asymptotically, because in this limit they com- 
mute with circulant matrices which all commute be- 
tween each other (using [M| and [3S]). A circulant 
matrix A with dimension n x n is defined as Ai-j = 



ai_jmodn- Therefore, we can introduce the notation 
k = (i — j) mod n which indicates the fcth diagonal 
of A. From [31] we state, that the eigenvalues of A are 
^Jrn = Efc=o Ofc e-*27rmfc/n^ m = 1, 2, ...n. If wc take the 
limit n — ^ oo the latter becomes the valid solution for the 
eigenvalue spectrum of all Toeplitz matrices (that belong 
to the Wiener class). As argued in Sec. Ill B the optimal 
input covariance matrix 7in is diagonalized in the same 
basis as the noise covariance matrix 7cnv Thus, 7in is 
asymptotically Toeplitz with quadrature spectra 



fe=0 

where jf^^^. are the fcth diagonal of 'y^^f (and the Fourier 
coefficient of a Fourier series) in the original basis. Since 
7in^(^) i^ Riemann integrable we conclude that 



1 



71^^ = ^/^^^ A; = 0,l,2,...,oo, (Bl) 

which provides the covariance matrix of the input state 
in the case of a global water-filling. 



Appendix C: Gauss-Markov process of order V 

Here we state the extension of the Gauss-Markov pro- 
cess to a Markov process of order V which is also called 
autoregressive (AR) process with white Gaussian noise. 
The underlying stochastic process is defined as [32] 

r 

= E(/.fcZt_fe + Wt, t = l,...,n (CI) 
fc=i 

where (jji, 4>2, (('v are the correlation parameters and 
Wt are identically and independently Gaussian dis- 
tributed random variables. This process is stationary 
(shift invariant) iff all roots of the characteristic polyno- 
mial 



p(y) = 1 - ^./.fe/ 



(C2) 



fe=i 



lie outside the unit circle \y\ = 1. If the process is station- 
ary, then the covariance matrix of the stochastic process 



(CI) is Toeplitz. Its spectrum is given by 



lARix) = 



Var(ZO 
1 - ELi (t>k e 



ikx 1 2 



(C3) 



where Var(Z'f) is the variance of Zt- Thus, if the noise 
correlations in both quadrature blocks in ( 26 ) are given 



by an AR process of order V one immediately can com- 



pute the capacity using ( C3 ) and the algorithm presented 
in Sec. HVj 
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